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Abstract 

We develop an analytical approach to scalar field theory on the fuzzy sphere 
based on considering a perturbative expansion of the kinetic term. This expan- 
sion allows us to integrate out the angular degrees of freedom in the hermitian 
matrices encoding the scalar field. The remaining model depends only on the 
eigenvalues of the matrices and corresponds to a multitrace hermitian matrix 
model. Such a model can be solved by standard techniques as e.g. the saddle- 
point approximation. We evaluate the perturbative expansion up to second order 
and present the one-cut solution of the saddle-point approximation in the large 
N limit. We apply our approach to a model which has been proposed as an ap- 
propriate regularization of scalar field theory on the plane within the framework 
of fuzzy geometry. 
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1. Introduction 



Roughly speaking, tiie term fuzzy geometry labels the noncommutative deformation of Rie- 
mannian manifolds which come with a Laplace operator with discrete spectrum. The idea is 
to truncate this spectrum and subsequently to deform the product on the truncated set of 
corresponding eigenfunctions in order to arrive at a closed algebra of functions. The most 
important example of a fuzzy space is certainly the fuzzy sphere, obtained by truncating 
the set of spherical harmonics Y/^ to the subset / < L together with the deformation of the 
product given by 

★ X-' - X-' * = [x"- * X^] ~ lEijkx'' , (1.1) 

where x*, i = 1, . . . , 3 are coordinates denoting a point on C R^. This deformation can 
be traced back, at least, to Berezin p]; the geometrical interpretation is essentially due to 
Madore p] . For a comprehensive review of fuzzy geometry, see jS] . 

The concrete physical motivation for studying fuzzy geometry is twofold: First, fuzzy 
spaces arise - similarly to other noncommutative spaces - in string theory rather naturally, 
when certain background fields are turned on [4J. Second, field theories on fuzzy spaces 
reduce in general to finite-dimensional matrix models and there are well-defined limits in 
which the classical field theories on the fuzzy sphere tend to the continuum field theories on 
the sphere and on both the commutative and the noncommutative plane. Therefore, fuzzy 
geometry might provide a powerful regularization scheme for quantum field theories, which 
would possibly have several advantages over a lattice regularization. 

A serious obstacle to using fuzzy geometry as a regulator is the fact that the taking of the 
limits which yield the commutative spaces does not commute with quantization: In [5], it was 
found that a finite "noncommutative anomaly" survives the limit from the fuzzy sphere to 
the commutative sphere. This anomaly in turn yields the well-known UV/IR mixing on the 
noncommutative plane [6] in a certain limit. An approximate analysis of the phase diagram 
of fuzzy scalar field theory |7] as well as corresponding numerical studies [8] confirm this 
point. However, modifications of the naive action used for quantization on the fuzzy sphere 
might reproduce the correct commutative limit [9]. 

To improve our understanding of the situation and to study the effects of the proposed 
modifications, an analytical handle on the model would be desirable. As mentioned above, 
field theories on fuzzy spaces are matrix models, and for these models, a large set of tech- 
niques has been developed over the last decades, which allows one, for example, to evaluate 
the partition function of certain matrix models exactly. Furthermore, gauge theory on the 
fuzzy sphere has recently been solved analytically by using localization techniques on a refor- 
mulation of the theory ^lOj. Unfortunately, scalar field theory turns out to be too complicated 
to allow for a complete treatment. 

Using the symmetry of the path integral measure, we can however show that the partition 
function of fuzzy scalar field theory is the same as that of a multitrace hermitian matrix 
model. Knowing this model explicitly would be very useful, as such a multitrace matrix 
model can be solved using e.g. the saddle-point approximation. As it turns out, a perturbative 
expansion of the kinetic term allows us to identify the multitraces appearing in the rewritten 
action order by order in a straightforward (but increasingly tedious) way. It is clear that 
such an expansion makes sense in a particular region of the parameter space of the theory. 
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Furthermore, this expansion corresponds to the high temperature (or hopping parameter) 
expansion which has been successfuhy used in the calculation of the phase diagram of scalar 
field theory on a two-dimensional lattice; for details, see [llj and references therein. 

In this paper, we develop the details of our technique and evaluate the perturbative 
expansion up to second order together with the solution of the thus obtained multitrace 
matrix model. Explicitly, we find the effective action 



^eff = 6 Y.^X^ - Xjf + ^(2,2) ( " A,)^) + ^4 ^^(A. - A,)^ + J] ^(A,) , (1.2) 



2 

I 

i>j i>j i>j i 

where Aj are the eigenvalues of the matrices encoding the real fields on the fuzzy sphere; 
the exact values of the coefficients '6(2,2) and .^4 can be read off from equations (|3.9p and 
()3.25p . The large behavior of the coefficients is given in (|4.4p and reads as 

^2 ~ — , ^2,2 = and ^4 ~ — • (1-3) 

Moreover, we apply the same technique to the proposed modification of fuzzy scalar field 
theory mentioned above, and compare the results to the ones from the undeformed model. 

The outline of this paper is as follows. We start with the construction of the model, com- 
ments on various special cases and symmetry considerations in section 2, before we present 
the actual perturbative expansion in section 3. The solution of the obtained multitrace ma- 
trix model in the large N limit is presented in section 4. In section 5, we consider the modified 
fuzzy scalar field theory and its corresponding multitrace matrix model, and we conclude in 
section 6. Our Lie algebra conventions, further useful formulae and some group theoretical 
background is found in the appendix. 

2. Scalar (/)^-theory on the fuzzy sphere 

2. 1 . The fuzzy sphere 

Consider the sphere S"^ = CP^ together with the Laplace operator obtained from the Fubini- 
Study metric. Its eigenfunctions are the spherical harmonics Yim with eigenvalues /(/-|-1) and 
degeneracy 2/ -|- 1. The quantized version of this space, which is known as the fuzzy sphere, 
is obtained by truncating the spectrum at a certain value L and deforming the product such 
that it closes on the corresponding truncated set of eigenfunctions. This yields a finite- 
dimensional algebra, which approximates the algebra of functions on 5^. To be more precise, 
note that in complex coordinates, all the with / < L can be written as homogeneous 
polynomials of degree L in both the homogeneous coordinates z^, a = 1,2, on <DP^ and their 
complex conjugates Za- As a basis for these polynomials, we choose the monomials 

^ai •••^«L% • (2-1) 

A suitable quantization of this basis is given bj0 

ai,...aijO){0\a(,,...ap^GnL(S)nl , (2.2) 



^This particular representation of functions on the fuzzy sphere can be found in |12) . 
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where |0) is the vacuum in the Fock space T of two harmonic oscillators with raising and 
lowering operators satisfying [ocajj] = baji- Evidently, Jii^ C T denotes the Hilbert space 
with L excitations and is its dual. The operators (|2.2|) take over the role of the spherical 
harmonics (j2.ip and form by construction a closed algebra. They are furthermore in a one- 
to-one correspondence with endomorphisms of sections of the line bundle 0{\)'^^ = 0{L) 
over CP^, and thus our quantization prescription is indeed the usual geometric quantization 
of CP^ 

Note that the operators (j2.2p act on the N = L + 1-dimensional irreducible representation 
of SU(2) and they can thus be translated into N x N matrices which become hermitian, if we 
consider operators corresponding to real functions. Thus, scalar field theories on the fuzzy 
sphere can be written in the form of hermitian matrix models. 

Furthermore, it can be shown that the Laplace operator A obtained from the Fubini-Study 
metric which is equivalent to CiCi with Ci = iSij^x^dk, where are Cartesian coordinates 
denoting points on 5^ C , turns into the second Casimir C2 in the representation obtained 
from the quantization procedure. 

Demanding invariance under the isometry group SU(2) suggests that after quantization, 
the integral over the sphere turns into a trace over the integrand. The appropriate normal- 
ization is given by the condition that dA = AttR'^. Therefore, 

IjAf ^ ^tr(/). (2.3) 

In the following, we will omit the hats over operators for convenience. 

Altogether, the fuzzy sphere is defined in terms of a matrix algebra approximating the 
algebra of functions on the commutative sphere and it is characterized by the parameters N 
and R. We will consider various limits of this matrix algebra in section 5.1. 

2.2. The model 

With the prescriptions of the last section, we can directly translate ordinary scalar (/>^-theory 
on the sphere to the fuzzy picture. Note that as in the commutative case, we scale the 
Laplacian by a factor of Our action is therefore given by [13] 

S = 7tr (-^c^,C2cI> + rcI>2+5(i)4) = ^tr (-^[L,,c^][L,,cI>]+rcI>2+5<I>4) , (2.4) 

where we introduced the shorthand notation 7 = . Furthermore, the Li are the genera- 
tors of SU(2) in the L + 1 = A^-dimensional representation of SU(2); see appendix C for our 
Lie algebra conventions. 

The partition function of the model reads as 

Za,r,g,N = / dfloi^) d = d/iD(^) Q ^ ^ " ^ ^ > ■ (2.5) 

The measure d^oi^) denotes the Dyson measure [T4] on the set of hermitian matrices of 
dimension N x N, 

d/iD(^) := '[[dRe{'^i,)'[[dlm{^ij) . (2.6) 

i<j i<j 



3 



This measure is invariant under the adjoint action of a unitary matrix 0, <I> — > Q^Q^; the 
normahzation is irrelevant for our purposes. However, we wih rely on the fact that under the 
decomposition 

$ = QAQ^ , A = diag(Ai,...,Aiv) , (2.7) 
the Dyson measure splits according to 

j diiom = /lldA, A2(A) / diJLHi^) . (2.8) 

i=l 

Here, d^ni^) is the Haar measure as defined in appendix A and A (A) is the Vandermonde 
determinant 

A(A) := det([Ari],,) = - A,) . (2.9) 

The explicit appearance of the radius R of the sphere is only interesting for approaching 
the planar commutative and noncommutative limits. We will therefore put = 1 in most of 
our subsequent discussion. 

2.3. The case a = 0; Hermitian matrix model 

If we put a = in (j2.4p . our model reduces to the well-known hermitian matrix model. 
This model is exactly solvable [15]. Not only can one straightforwardly reduce this model 
to a model of eigenvalues, thereby reducing A^^ integrations to A^, but one can even give an 
explicit expression for the partition function in terms of orthogonal polynomials. 

Using the decomposition <I> = OAil^, the dependence on Q. drops out and the model turns 
into 

N 



Za=, = I HdA, A2(A) I d^,Hm e-^^'^(^^^'^+^^'' 
1=1 
r ^ 



1=1 

which is the description of a one-dimensional gas of eigenvalues in a quartic potential with a 
two-body repulsion. To obtain the exact partition function, consider the monic polynomials 
Pn of degree n normalized according to 

I dAp^p„e-^^*'-(^^'+3^') = h^6mn . (2.11) 

Then the partition function reads as 

Z = N\ Yl hi . (2.12) 

i=0 

As it is well-known, the large A^ limit of this model has a third order phase transition at the 
curve 

7 2 2 1 Q^ 

q = r = — . 2.13 

for negative values of r. At this curve, the double well potential becomes sufficiently deep, 
so that the eigenvalue distribution splits into two components. This phase transition will be 
expected to appear also for a / 0, gued in [7J. 
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2.4- The case a ^ 0: The full fuzzy field theory 

With a 7^ 0, our model ()2.4p becomes considerably more difficult than the hermitian matrix 
model, as the "external matrices" Lj do not commute with unitary matrices, and thus prevent 
us from performing the integration over the angular degrees of freedom encoded in Q. 

In [16] , the partition function for the hermitian matrix model with an additional external 
matrix A and the interaction term V{A^) has been calculated exactly by counting graphs; 
this solution was reproduced in [17J using the character expansion method. Unfortunately, 
having three external matrices Lj as in our model ()2.4|) renders this method essentially useless. 

However, one can make the following statements: The kinetic term does not depend on 
tr ($) since C2I = 0. Therefore, when J d;^i^ tr ($(72^*)) is rewritten as a function 
of the eigenvalues of it only depends on the differences (Aj — Xj). Since there has to 
be a permutation symmetry between the different eigenvalues as well as invariance under 
Aj — > — Aj, the kinetic term is in fact a function of the form 



Besides yielding a nice consistency check for our calculations later on, this form offers a 
physical interpretation: The kinetic term corresponds to an additional attractiv^ force bet- 
ween the eigenvalues, which will compete against the repulsive force stemming from the 
Vandermonde determinant. 

Altogether, we can make the following qualitative predictions: We expect three different 
phases, each characterized by the dominance of one of the parameters a, r and g. 

First, there clearly has to be a remnant of the disordered phase of the matrix model 
with a = 0, in which g is dominant. In this phase, (| tr {^)\) is small due to the eigenvalue 
repulsion. As the eigenvalues are confined in a single well potential around the origin, also 
(tr ($^)) is smah. 

Second, there can be a non-uniform order phase, in which r strongly influences the shape 
of the potential and the double well becomes sufficiently deep. Here, (| tr ($)|) is still small, 
but (tr(<I>^)) becomes considerably larger. 

Third, there should be a uniform ordered phase, in which the kinetic term plays a signifi- 
cant role. As the force introduced by this term corresponds to a potential which vanishes for 
$ ~ To, we expect that in the large A'' limit, the expectation value of (| tr ($)|)^ ~ (tr ($)^). 
Surprisingly, such a phase seems to be present even in the absence of the kinetic term [18] . 

2.5. From fuzzy scalar (j)^ -theory to multitrace matrix models 

It is easy to see that any hermitian matrix model can be rewritten in terms of multitrace 
expressions: Since the Dyson measure is invariant under the adjoint action by an arbitrary 

^The force cannot be repulsive, as the scalar field theory on the fuzzy sphere with g — Q and r > is 
evidently stable for finite A''; it can, however, be repulsive on short distances. 




where 





5 



unitary matrix 17, 

d/UD(^) = dfioi^'^^^) , (2.16) 

the part of the action relevant to integration has to be invariant under this action, as weh. 
Formally, we can write: 

J df^om e-^ = J df^om e-^=ff with e-^^«[*l = ^^^^^^ J df^ni^) e-^[^*^'l , 

and Scs is clearly invariant under $ — > ri^ri'''. Now the only functions of ^> invariant under 
this transformation are products of traces, and we thus have 

SeS = 5]5„tr(cl>")+^>.„„tr(cl>")tr($")+ ^ w tr tr ($-) tr (c^.'^) + . . . . 

n n,m n,m,k 

Such an action is certainly much easier to study than the original one involving an external 
matrix. Our goal will therefore be to turn our model of fuzzy scalar (/>^-theory at least 
approximately into such a multitrace matrix model. 

Multitrace matrix models became popular at the beginning of the 1990ies, as they were 
found to arise naturally when higher order curvature terms were included in matrix mod- 
els describing two-dimensional gravity [19], see also [20j. Including these curvature terms 
amounts to modifying the kinetic term of the matrix model by adding tr (A^A^), where A 
is a fixed external matrix. Before this model was solved exactly in fTB], the authors of [19j 
approximated it by the simplest multitrace matrix model of the form 

S = atj:{^'^)tr{^'^)+rtv{<^>^)+gtT:{^^) . (2.17) 

This matrix model has been considered as a model for string theories with c > 1 and ad- 
ditional curvature terms like / d^a^E?. It was found that the phase diagrams of both 
theories share many features |20j . 



2.6. The most general hermitian matrix model 

As a remark, note that one can push the analysis of the structure of hermitian matrix models 
using the symmetries of the measure even further. The full symmetry group of the Dyson 
measure d/i£)($) on the ensemble of N x N hermitian matrices is in fact SO(A^^), as 

duDm /(CD) ~ d^'^^' /(cD%) . (2.18) 

where d^^^'^ is the ordinary volume element on . This implies that the only relevant 
part of the action S is the one which is invariant under SO(A^^) rotations of the components 
5'eff) and therefore, the action has to be a function of tr ($^): 

fdfiDme-^ = /"d;UB($) e-^^« with 5eff = s„ (tr ($2))" . (2.19) 

In this expansion, however, the large N limit becomes entangled with the rewriting of the 
model in terms of tr (^D^), which renders the reformulation useless for most purposes. 
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2. 7. The toy model: Hermitian 2x2 matrices 

For simplicity, let us consider for a moment the case L = 1. The matrices $ act here in the 
fundamental representation of SU(2), i.e. $ is a hermitian 2x2 matrix. We expand <I> = ^^'^r^ 
in the generators = (tq, t^) of SU(2); see appendix C for more details on our conventions. 
Together with the commutation relations [Ta,r;,] = ^/2i£abcTc, the action (j2.4p simplifies to 

S = 7(2a$''$'' + r$''$'' + 5tr(($^r^)^)) . (2.20) 

Since = tr ((<!> — ^1 tr (<I>))^) = tr (<1)^) — ^ tr (<!)) tr (<I>), we can rewrite this expression 

in terms of traces: 

S = 7((2a + r)tr(^>2) -atr(^>)tr($) +c/tr($'^)) . (2.21) 

As in the case of the pure matrix model, one can decompose $ into the product rJAfi^, 
where is a unitary matrix and A = diag(Ai,A2). The Dyson measure d/i£)(<I>) will again 
split into the measure on the space of eigenvalues dAidA2 and the Vandermonde determinant 
(Ai — A2)^. In total, the partition function for L = 1 reads as 

Za,r,g = J dAidAs (Ai - Aa)' ^-Piiaix^-x,)^+rixl+xl)+g(xi+xl)) (2.22) 
We will use this special case for comparison with our perturbative results for general N. 

3. Perturbative expansion of the model 

3.1. Motivation 

As mentioned in the previous section, our goal will be to turn fuzzy scalar field theory into 
a multitrace matrix model. Since we do not know of any direct route to integrating out the 
angular degrees of freedom exactly, we will have to resort to perturbative techniques. This 
approach is also suggested by the lattice analysis of scalar (/)^-theory on , where a hopping 
parameter - or high-temperature - expansion has been successfully used, see e.g. [H] and 
references therein. 

Assuming that the coefficient of the kinetic term, a, is small - or equivalently that the 
coefficients r and g are large, it clearly makes sense to expand the model in powers of a. To 
demonstrate this, we use the high-temperature expansion to compute the specific heat of our 
model for N = 2. The results are shown in figured! where the specific heat ^ is plotted in 
the case N = 2 for both the pure matrix model and the full theory with high temperature 
expansion up to 0{a^). In the pure matrix model, one clearly sees that the critical line of the 
third order phase transition is a parabola. The plot corresponding to the full model shows 
that the critical line for the second order phase transition is indeed a straight line, and the 
slope extracted from this data deviates from the numerical prediction by only 3.25%. 

Note that the matrix model phase transition is certainly also present in the full model, 
however, it is invisible in the second plot, as the scale of the maximum is far larger than 
the scale of the parabolic valley. In general, for large g, the specific heat 'rf := -^logZ 
interpolates between ^ at large positive and negative r. However, there is an intermediate 
phase, where ^ approximates j which corresponds to the negative r phase of the pure matrix 
model, see figure [2j 
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Figure 1: On the left, the specific heat := ^ log Z for the pure matrix model with o = 
is plotted. The right plot presents the result of the 0(a^)-approximation to fuzzy (/)^-theory 
for = 2. 
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Figure 2: Cross-section through the right plot of figure 1 at 5 = 50. 
3.2. Perturbative expansion - 0{a^) 

We will use again the Gell-Mann basis and write $ = $'*r^. Note that the kinetic term only 
depends on as C2 annihilates tq ~ 1. We start from the expression 

g/37a3.-X„,$'' ^ 1 + pja^^Kab^'' + ^^^^"^ $°Ka&^^ ^"K^d^'^ + ... , (3.1) 

where 

= tT:{[Li,ra][Li,n]) . (3.2) 

As in the case of the pure matrix model, we diagonalize $ = ^}AD) using the unitary matrix 
Q. This matrix contains the angular information on $ and has to be integrated over at the 
end. We thus have 

= tr (r^QAfit) , (3.3) 

and we need to compute 

J d/Xi/(J])i('a6tr(T"0A0t)tr(r^QA0t) (3.4) 
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at first order and 

j diJLH{^)KabKcd tr {T^nKQ^) tr (r''J7 AO^) tr (r^fiAQ^) tr {r'^nKQ^) (3.5) 
at second order. 

After rewriting these expressions as sums over traces in irreducible representations, we 
can apply the orthogonality relatioijf] 

j d^,Hm [pmh [p\^)]ki = -^[^^^ii^jk , (3.6) 

to perform the angular integration. We therefore continue by rewriting ()3.4p as 

ii {{T^nm^) ® {T^nm'^)) = tr((r''®r'')(l^®17)(A«)A)(f]t0f^t)) =. tr([a8)6]). 
The tensor produc10 (8) [i] splits into \^ © ^, and thus we have 

tr([a®6]) = tr03([a®6])+ trra([a®6]) . (3.7) 



The integral is now easily performed by applying the orthogonality relation (j3.6p . and we 
arrive at 



(3.8) 



Note that tr p(A (8> A) is actually the character of A in the irreducible representation p. The 
traces can be easily computed and the results are given in appendix B. With these formula, 
we obtain 



j dnHimab^''^' = trK(^^^(trA)2 + ^^trA2) 

= ^A^^(trA^-^(trA)^) = iiV^^.-A,)^ = 



(3.9) 



where we used 

L 

ii{K) = tr(C2) = Y,^2l+l){l{l+l)) = \L{l+Lf{2+L) = \{N-1)N'^{N+1) . (3.10) 
and 

E^'-fE^^^^- = (I-f)E^'-fE^^^^- = ^Y.^>^1. - '^hX, + \]) . (3.11) 

i i,j i i>j i>j 

One easily checks that (13. 9p is indeed independent of tq: for A = l^r and, equivalently, 
for eigenvalues Ai = ... = Xn, the expression vanishes. Note that we can replace A by $ 
everywhere, keeping the expression independent of the angles Q. Therefore, one can "invert" 
the integration and replace Kab^"^^^ in (j3.ip under the integral directly according to 

j dpDmKab^"^' = JdfiDi'^)N^(^tri<^>^)-^tr{^f^ . (3.12) 
^For a proof of this relation, see appendix A. 

^Note that our notation here obviously deviates from the standard one. For more details on this, see 
appendix B. 
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3.3. Perturbative expansion - O(a^) 

The calculations for higher orders in a and thus K are more and more involved, but never- 
theless straightforward. For our approximations, we would like to have one more order in a, 
so let us calculate the involved terms: 

J dflD{^)Kab^''^''K,d^^^'^ , (3.13) 

which turns into 

j d/iD($)X]^j^^^^<Jfe^cdtrp(r°0r''0r=®r'^)trp($O$0$0#) . (3.14) 



after using again the orthogonality relation (13. 6p . 

We start by considering the terms KabKcd p{t"' ® t'' r'^ (8> t'^). The traces of this form 
break up nicely into ordinary traces and multitraces as above, and the detailed formula are 
found in appendix B. Additionally, we use the identities 



if,,K,rftr(rVW) = KabK,4l^6''''6'U^d-V- 



(3.15) 



which are easily proven with the identities given in appendix C. 

We have now the following contributions in the various representations, multiplicities 
already taken into account: 

^ ■ ^ ((^ + I) + + ^) + d^"' + d-^d'"^^) KabK,^ 

^ • "K(^ ^ ^) ^ ^ ^) ^ [kd-'^d'^''' + d^^'d'^-) KabK,^ 

o 

-\{^[l-^^{iTKf+{2- 1^ trET^ _ (^^^abe^cde^^ace^bde^^ K^bK,^ 
^ ( (l - ^) ( tr K)2 + (^2 - 1) tr _ (^l^abe^cde ^ ^ace^bde^^ K^bK,^ 

(3.16) 

Note that - as expected - the sum over the above terms vanishes. To obtain the correction 
to second order in a, we have to divide each term by the dimension of the corresponding 
representation and multiply by the character of A (or <I>) in this representation. The latter 
are taken from table ()B.8|) . and putting all together, we obtain 

j dllDrnKab^^^^Kcd^'^'^ = 

f (2tr + itrKf)iti + Stitg - 2(tft2 + U)N + tjiN^ - 6) , 

-5tf - Atit3 + 3tl + (lOtjta + U)N - 2{2tih + tl)N'^ + t^N^ 
^ iV(-36 + iV2(_7 + ^2)2) 
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where we used the shorthand notations 

K^K := (^^d'^^^d^^^ + d^'^d'''^^^ KabK,d (3.18) 

and 

:= {t,y := (tr(A^))^ (3.19) 

As discussed in section 2.4, we can in fact rewrite the integrand in (j3.17p in terms of the two 
functions: ^ 

H4 = Y.{Xi-\j)^ and El = (^(Ai-A,)2) . (3.20) 

i>j i>j 

For this, note that 

2 ^(Ai - XjY = Ntg - Q tg-ltl + • • • - ^ ^) htg^l + Ntg 

i=0 



i>j 



(3.21) 



In our case, iVt4 uniquely identifies the contribution of H4, while N'^tl marks the contribution 
of H| . We thus read off that 

where 

ai = -2, 02 = 1 , /3i = (7V-2)(iV-3) + 5(iV- 1) , = -5 ; (3.23) 

performing the complete calculation as a consistency check yields the same result. Using 
additionally that 

(tri^)2 = {trC2f = l{N-lfNHN + l)\ 
trK^ = (trC|) = l{N-lfN^{N + lf, 
we can further simplify our result to 

I df^nimat^'^'^'K^d'^''^' = I d^^D{<^) (-2S4 + H^)+ ^^^^^ 

+ Ar(-36 + iV2(_7 + ^r2)2)(/3iS4 + /32Hi)K-i^ . 

5.^. Formula to all orders 

In general, every order in a of the expansion of the kinetic term can be rewritten as a 
polynomial of the eigenvalues in this manner. Multiplying by the Vandermonde determinant 
only increases the degree of the polynomial. 
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To all orders, the expansion reads as 

i n=0 pey(2n) ^'^'^ (,3.2bj 

X trp(r"i r^i © . . . (g) r"" (g) A)xp(A) , 

where y(2n) is the set of irreducible representations of SU(A^) corresponding to Young di- 
agrams with 2n boxes, A is the diagonal matrix of the eigenvalues of ^ and Xp(A) denotes 
the formula for the character in the representation p, evaluated at A. 

If one were to use orthogonal polynomials for solving the approximated model, note that 
the only integrals which would appear are of the form 

Mia / nr /i_L^\ /i_L^i^g2^ 



3 + n \ / 3 + n 3 
4 4 '2'4r 

and since they are thus exactly soluble for ^ > 0, which corresponds to our case, one can 
in principle perform exact analytic computations at every order in a. In practice, the terms 
appearing will be more and more involved and it is thus advisable to use a computer algebra 
program. 



3.5. Verification for N = 2 



To verify the above description, let us briefly compare to the exact results in the case N = 2, 
i.e. L = 1. Prom (j2.2ip . we obtain 

J dfiDi^) ^""Kab^' = J dfiDi^) 2tr($2) _ tr($)tr($) , (3.27) 

which is identical to ()3.12p in the case N = 2. 

In the case of the term at second order, the general expansion of the tensor product is 
truncated according to 
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(3.28) 



As a consistency check, one can sum tr^ over these representations which yields the usual 
tr (r^) tr (r^) tr (r'') tr (r'^). Using the identities (vahd for L = 1) 

i^,fci^e,tr(TVW) = \b''H'''KabK,^ = U^vKf, 
and the formute given in table ()B.6p . one easily verifies that 



(3.29) 



j dllHi^) $"i^,b$^$=A'erf$^ = j d^H(^) {-At^^^ + Q{tT^^f -{tl^Y 

in agreement with the term of order 0{a?') in the Taylor expansion of ()2.22p . 



(3.30) 
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4. Solution of the model for large N 



4-1. Large N limit 

Numerical results [8j suggest, that the parameters of fuzzy scalar field theories at different L 
can be rescaled to make the corresponding phase diagrams collapse to one unique diagram. 
This implies in particular, that we can study the large limit, in which the formulas simplify 
and we can apply the saddle-point approximation. Let us therefore extract the large N limit 
in the following. 



Using the identity r^^r^; = Su^jk and 



(4.1) 



(4.2) 



Ka,K,d tr {t'^tW) = \ tr ([L„ r<^] [L„ r*]) tr ([L„ [L,,A) tr {t'^tW) 

= itr([L„r'^][L„rn)tr([L„r''][L„rn)tr(r^r-r''r-) , 

we find that 

i^afe^cdtr(rVVr'^) = \nHi {UUL.Lj) + \N{iT {ULi))\ 

We have furthermore LjLj = ^(iV^ — l)lAr and tr(LjLj) = ^N{N'^ — l)5ij, from which it 
follows that 

^NHiiULiLjLj) + ^N{tTiLiLi)f = lN'^{N^-lf, 

itr([L„L,][L,,L,]) = -N{N'-1) . 

Thus, the dominant term in (j3.16p is ( tr K)^ ~ N^, in particular the term K^K is negligible 
in the large N limit. Altogether, (jS.lTp eventually reduces to 

j d^iDimab^^^'K^a^'^" = y"dMD(^)(-^H4 + ^H2) , (4.4) 

where the second term is precisely the next term in the Taylor expansion of the exponentiated 
correction in 0{a^). The first term is new, and corrects the matrix model. Note that it has 
the right sign to be exponentiated again and keep the partition function well-defined. 

One is clearly tempted to conjecture that at reth order in the high-temperature expansion, 
the term X]j>j(Aj — Aj)^" is new, while all the other terms just correspond to terms in the 
Taylor expansion of corrections at smaller orders. This, however, seems not yet possible to 
prove with the techniques at hand. 

If this assumption is correct, we can reexponentiate the terms and arrive at the following 
matrix model, written in terms of eigenvalues: 

i i>j 

Note that even if the assumption was not correct, this would still be a useful approximation 
valid up to order C?(a^). 

In the large N limit, we rescale Xi ^(t^) — A (a;) with < x < 1 and the sums are 
turned into integrals: X^^^ N dx. We thus arrive at the action 

S = jN [ dx(rX\x)+gX^{x) + N I dy{ - ^N{X{x) - X{y)f 



+ ^4n\X{x) - X{y)f - i In |A(x) - X{y)\)) . 
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The saddle-point approximation 

For both single and multitrace hermitian matrix models, the saddle-point approximation 
provides a powerful tool for analyzing the phase structure. However, one has to integrate 
out the angular variables first and reduce the matrix model to the eigenvalue formulation. 
This is due to the fact that the angular degrees of freedom are zero modes in the partition 
function and keeping them would render the saddle-point approximation invalid. 

In the large limit, we can define scaling limits for the coupling constants, such that all 
of the terms contribute. From the log-term, we see that we have to pull out a factor of A^^. 
Supporting the validity of our approach is the fact that we can choose a scaling which agrees 
with the scaling found numerically [8j: 

a = N^^a , r = N^'-f , g = N^^g and X{x) = iV^^A(x) , (4.7) 

where 

Oa = -k, Or = I, Og = I and = 3 . (4.8) 
Note that the correct conversion is the following: 

Ob = -i + er-ea + 2ex and = -i + Og-iOa + ^ex , (4.9) 

where h and c are the coupling constant in the normalization of [S]. 
Altogether, we have the following partition function: 

&X exp(-iV25) , (4.10) 

where 

"1 . rl 



S = 47T dx(fX\x)+fX\x)+ / dy(-f(A(x)-A(y))2 

Jo ^ Jo ^ ^ (4.11) 



_|_ 47Va? '~ 



(A(x)-A(y))^-^ln|A(x)-A(y)|)) . 



With our scaling, the saddle-point method becomes exact in the large limit. 

Rewriting the integral over y in (j4.1ip as an integral over X{y) introduces the eigenvalue 
density u(A) = dy/dX, and the moments 

cq = J djl u{jl) = 1 and '^^ ~ j ^(A)/^* ' ^ ~ (4-12) 

appear. Since our problem is symmetric with respect to A ^ —A, it is natural to assume that 
the eigenvalue density will possess this symmetrjjf] and therefore the odd moments ci and C3 
have to vanish. Together with suppressing constant terms in A, this reduces the action (j4.1ip 
to the following one: 

S = 47r^'dx(^(f-f + 6fic2)P(x) + (5 + ^)A4(x)-^^'dyln|A(x)-A(y)|) , 

(4.13) 

which corresponds to the action of the minimal multitrace matrix model studied in with 
shifted coefficients. 

^The results of [18] indicate that this is not necessarily so; however, we will postpone the detailed discussion 
of this point to a future publication. 
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4-3. The solution 

Since the solution to our model can be derived by well-known methods, let us be brief in the 
calculation. Assuming that we have a one-cut solution, for which u{X) has support on the 
interval [—5,5], the saddle-point equation, which is the equation for the stationary points of 
(|4.13p rewritten using the eigenvalue density u{X) under the constraint cq = 1, reads as 

47r - f + ^ca) A + Svr (5 + ^) = j^^ d/i , (4.14) 

where the integral on the right hand side is a principal value integral. The solution to (j4.14p . 
obtained in the usual manner, is 



n(A) = [Ar-a + 12™' ca + A (g + ^) 6' + 8 [g + ^) X') ^ 6'^ - X"^ ; (4.15) 



the familiar polynomial deformation of the Wigner semicircle distribution. The normalization 
condition cq = 1 yields 

^ttS"^ {4r + ISgd"^ - a + S-Ko^ {4c2 + 36'^)) = 1. (4.16) 

Solving this equation for C2 results in 

2 + 7T6^{d-Af-9{2g + d\)S^) 
= 12^^W ' ^^-^^^ 

and using the original definition of C2 yields the self-consistency condition 

TT6'^i8f + 36g6'^+d,{3d,7r6^{8 + ir{2g + d'^TT)5'^)-2))-4 = 0. (4.18) 

4-4- The phase diagram 

As discussed before, we expect a non-uniform order phase in the phase diagram, and therefore 
a region in the parameter space for which the one-cut assumption is no longer valid. In this 
region, the eigenvalue density (j4.15p will become negative and therefore ill-defined. Looking 
at the polynomial part q{X) of u(A), we see that the transition should occur, when q{0) = 0. 
Together with the normalization of cq and the self-consistency condition for C2, this gives a 
critical curve C marking the critical values of g depending on the other parameters. Explicitly, 
C is defined by the following two branches: 

= ^ (^-63a2 - 8af + 16r^ ± {4r - d)y/l6r'^ - 8df - OSa^) . (4.19) 

The same result can be obtained straightforwardly in a different manner: Let us again 
assume a symmetric eigenvalue density which allows us to turn to the action 

S = S2tr{<^>^)+S2,2itr{<^^)f +S4tr{<^>^) . (4.20) 

Furthermore, we define the effective coupling S2 = S2 + S2,2(ti' (^^))- We can now treat our 
model as if it was a pure, single-trace hermitian matrix model with §2 as the coefficient of the 
quadratic term. In the regime in which the eigenvalue density has support on two disjoint 
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cuts, we have the usual formula (tr (<1>^)) = Using the equation §2 = —2^JNsi for the 

one-cut two-cut transition, we are lead to the rescaling of fields we used above and eventually 
arrive at the curve C as a self-consistent solution. 

To compare this phase transition to the ones obtained numerically, recall that in [8], 
three distinct phases were found. The transition between the first and the second phase is 
of thircfl order, i.e. the specific heat has a discontinuous first derivative at the curve roughly 
described by equation (j2.13p . The transition between the second and the third phase is of 
second order, i.e. the specific heat diverges, and numerical results suggest that the critical 
curve is approximately given by the equation 

g = (0.19±0.08)\/]Vr . (4.21) 

The phase boundary we find is a deformation of the matrix model boundary, cf. fig. 
[31 In this deformation, the region in which the one-cut solution is valid is increased by 
turning on a. This is in agreement with the predictions derived from the numerical results, 
that the triple point at which the third-order phase transition meets the critical line of the 
second order transition is located at larger values of r than in the case of the pure matrix 
model transition. Therefore, our results suggest to interpret the turning point of the curve 
given by as the rough location of the triple point, and we will use this interpretation in 
the following section. A comparison with the numerical values confirms this interpretation. 
The turning point is located at {r,g) « (—2.7,0.25) which is to be compared to the result 
(r, g) = (—2.3 lb 0.2, 0.52 it 0.02) from the numerical computations [8j, a reasonable agreement 
for a second order approximation. 

We expect the critical curve (|4.2ip to emerge in a future discussion of the full solution 
space to our multitrace matrix model. 




2468 10 2468 10 



Figure 3: The one-cut-two-cut transition for the matrix model and the boundary for the 
region of validity of the one-cut solution to the multitrace matrix model for a = 1. 



A detailed calculation proving this will be presented in a future publication. 
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5. Modification of the model 



5.1. N ^ CO limits of the model 



The various limits with N ^ oo for scalar field theory on the fuzzy sphere were discussed 
in [5]. To make these limits more transparent, note that we can obtain a noncommutative 



coordinate on Sp using the noncommutative analogue of the coordinate x* ~ z^a^^gZ^: 



= ai^^M=aR with x^x^L = -R^IL , (5.1) 

^iV2 _ I ' Wn-1 Wn~i ' ^ ^ 

where TCn-i is the two-oscillator Hilbert space with L = N — 1 excitations. These coordinates 
satisfy the algebra 

[x\x^] = i-j===eijkx'' . (5.2) 
There are now the following large N limits: 

> The limit of the commutative sphere 5^, for which we simply take — > oo with fixed 
radius R. 

> The noncommutative planar limit, in which we consider the tangent plan^ at the north 
pole. That is, x^ — > i? and therefore the correct limit is iV^ — > oo and R"^ = which 
results in the noncommutative coordinate algebra [x^,x^] = i6 on Rg. 

> The ordinary planar limit, for which N ^ oo and R^ ~ A^^"^ with > e > 1. In this 
limit, the coordinates x^,x^ turn into the commutative coordinates on IR^. 



5.2. Commutative and noncommutative (j)^ -theory on the plane 

As mentioned in the introduction, the naive action for scalar field theory on the fuzzy sphere 
will not reproduce standard (p^-theovy on R^. There is a UV/IR mixing [21] even in the case 
of the fuzzy sphere, which is, however, finite for every finite value of L [5]- In [9], it was 
shown that the deviation from the commutative theory is given by tadpole diagrams, which 
can be removed by modifying the action in the following way: 

S = 7tr (-^$(C22L(C2) + t-|i2(L,t))$ + rcl>2 + 5^^) , (5-3) 

where Zl[C2) is a power series in C2 corresponding to wave function renormalization, and 
R{L,t) is a function logarithmically diverging for large L. 

The numerically obtained structure of the phase diagram of fuzzy scalar field theory [8] 
suggests a comparison with that of (/)^-theory on , as the phase transition corresponding to 
the line is also seen there: The existence of a second order phase transition in scalar c/i^-theory 
on 1R2 has been rigorously established in [22] . This phase transition is due to spontaneously 
broken refiection symmetry (p —cj). The exact structure of the phase diagram was calculated 
only later using lattice and other techniques [23j . It was found that for the model 

Ce = \{V(t^? + \m'cf' + \\<l)^ , (5.4) 
more detailed discussion would use complex coordinates and the stereographic projection. 
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the phase transition on a lattice occurs along the line A = — 10.24(3)m^ in the A-m^-plane. 

Note that the third order phase transition stemming from the matrix model part together 
with its triple point with the second order phase transition is absent. This phase finds, 
however, a nice interpretation in terms of i;^^-theory on the noncommutative plane. There, 
the theory exhibits a new striped phase, which breaks translational invariance. This phase 
has been predicted in [24] for d > 2 using a self-consistent Hartree treatment and in [25] via 
a one-loop renormalization group analysis. Numerical evidence for the existence of such a 
new phase has been found in |26j . As argued in [27], appendix B, a definitive confirmation 
of this result is still pending. 

Altogether, we conclude that naive (/)^-theory on the fuzzy sphere serves most likely as 
a regulator of (/)^-theory on the noncommutative plane. Furthermore, we expect that in the 
modified model ()5.3p . the triple point moves off towards r = — oo, which would be a clear hint 
for a fuzzy field theory tending indeed to quantized (p'^-theory on the plane in the appropriate 
large N and R limits. 

5.3. Modified fuzzy (p'^ -theory 

Let us now study the effect of the wave function regularization on the phase diagram, as- 
suming that Zl{C2) ~ 1 + kC2. That is, we consider the modified action 

S = 7tr (-^$(C72 + KC2C2)$ + r$2+^$4^ . (5.5) 

This modification is particularly simple to realize in our formalism, as it amounts to replacing 
the kinetic matrix 

Kab Kab ■■= Kab + uKacKcb ■ (5.6) 

We thus have the following new quantities: 

tr (K) = liN - 1)N\N + 1)(3 + 2{N^ - 1)k) , 

tr(K2) = ±n'^{N^ -if {10 + 5{3N^ -4) n + 2{3N^ -9N^ + 8) K^) . 

To have both contributions in K survive the large limit, k = N^"-k has obviously to 
scale with 0^ = —2. After rescaling, we have K, and one can easily show that K^K is still 
negligible compared to the contribution from (trAT)^. A straightforward calculation shows 
that a is simply replaced by 

a = a{l + |k(a2 - 1)) (5.8) 
in all equations, and we therefore have in the large N limit 

~a = a{l + |k) . (5.9) 

As discussed above, increasing a moves the triple point away from the origin, so we verified 
that introducing k has the desired effect. 

6. Conclusion 

We gave an explicit algorithm for turning the kinetic term of a fuzzy scalar field theory 
into a multitrace expression which behaves equivalently under the functional integral. This 
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algorithm uses a perturbative expansion of the kinetic term, and it produces the exact result 
in a - in principle - straightforward manner at every order in the perturbative series. 

Explicitly, we evaluated the multitrace terms for scalar i;^^-theory on the fuzzy sphere 
up to second order and presented the one-cut solution in the large N limit. We found that 
this solution gives rise to a consistent deformation of the region of validity of the one-cut 
assumption in the case of the pure matrix model. Furthermore, this deformation exhibits 
a turning point, whose location is in good agreement with the location of a triple point 
conjectured from the numerical data. 

Eventually, we studied a modification of the original model which is believed to provide the 
correct regularization of two-dimensional ^''-theory within the framework of fuzzy geometry. 
The multitrace model arising from the deformed theory is a deformation of the multitrace 
model previously obtained and the observed changes in the phase diagram suggest that the 
modification is indeed useful. 

In a future paper, we will study the full solution space to the multitrace matrix model for 
both the naive and the modified 0^-actions. Furthermore, we will try to apply the technique 
of orthogonal polynomials for solving the model exactly also at finite N. Also, we want to 
probe the multitrace model numerically in order to gain solid data for comparing it with the 
full fuzzy field theory. 

Further directions for future research are also evident: As our rewriting of the kinetic term 
is independent of the potential, one can easily solve further models with different potentials, 
as e.g. </)^-theory or even sine-Gordon theory on the fuzzy sphere. Also, field theories on 
more general fuzzy spaces such as higher-dimensional complex projective spaces [28], fuzzy 
flag manifolds [29] and fuzzy projective algebraic varieties [30j are now accessible analytically. 
In particular, one should study which aspects of non-renormalizability of e.g. (/)^-theory in 
four dimensions remain, when the theory is considered on fuzzy <DP^ x <DP^ or fuzzy CP^. 

After scalar field theory and Yang- Mills theory on the fuzzy sphere have been approached 
analytically, it might also be interesting to consider the matrix model corresponding to general 
relativity on the fuzzy sphere, which has been constructed in |31j . 

Eventually, let us state a rather bold conjecture: Multitrace matrix models were first put 
forward and studied in the context of trying to define string theories with central charge 
c > 1 [191 [20] . It is well-known that such theories exhibit a phase transition with complex 
critical exponent, an earmark of high instability, and this instability is believed to correspond 
to a crumpling of the Riemann surfaces. There is a natural class of multitrace matrix models 
which arises from applying our approach to scalar field theories on the fuzzy sphere and fuzzy 
Riemann surface^ of higher genera. Furthermore, the truncation of the algebra of functions 
provides a rather natural mechanism for keeping the Riemann surfaces from crumpling in an 
uncontrolled way. It might indeed be worthwhile to look for a connection between fuzzy field 
theories and the long-sought definition of c > 1 string theories. 

*For a definition of such spaces, see e.g. [32]. Alternatively, one can define them by considering an embed- 
ding of the Riemann surfaces into CP" and using the techniques of [30] . 
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Appendix 

A. SU(A'") orthogonality relation 

In the above discussion, we used the orthogonality relation 

J dfiHi^) [pmij [pH^)]ki = , (A.i) 

where Q G SU(A^), p is a finite-dimensional, unitary, irreducible representation and denotes 
its complex conjugate. The measure dyU_H-(0) is the Haar measur^ on SU(A) normalized 
according to / dfiHi^) = 1- 

We briefly recall the proof of this relation, as given e.g. in [33]. Consider the integral 
I = f dpni^) p{^) Ap^ {Q) . Using the identity p{T)Ip^{T) = I, which is a trivial consequence 
of the Haar measure being invariant under translations — > TO, T G SU(A), one can 
show that p(T)/ = Ip{T). Thus, we conclude by Schur's lemma that I = cl. Choosing 
Amn = Smjdnk, wc obtain / dpH{^)[pi^)]ij[pH^)]ki = Cjkdii- Tracing this expression over 
the indices i,l yields that dim{p)cjk = 6jk, from which (jA.ip follows. 



B. Trace formults 

In this appendix, some explicit formula for traces of Kronecker products of matrices 
r'' . . . projected on certain irreducible representations p are given. Such traces will be 
denoted by trp. Furthermore, we use the shorthand notation 

(^1 . . . fimKi^i ...un)... := tr (r^^ . . . r'^™) tr {t^^ ...t^-)... . (B.2) 

Curly traces {pa . . .} denote a sum over all distinct permutations of the enclosed factors. 
Note that simplifications arise for traceless factors. It should be stressed that our labeled 
Young diagrams are not the usual ones: The labels in the boxes in the Young diagrams 
describe the action of the Gell-Mann matrices corresponding to the labels. For example, 
means the action of ® in the representation □□, which reads as 

'^^ijM ■■= \{jik'^ji + r'^k<i+^ii^^k + ^ji<kj ■ (B.3) 
®i.e. the measure which is invariant under left or right muhiplication by /5(T), T £ SU(A'^) 
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The tensor products of the fundamental representations we are interested in decompose in 
the following way into irreducible representations: 
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(B.4) 
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or, for convenience also written in terms of dimensions: 

N'^ + N N'^ -N 
N -N = h 



N-N-N -N 



2 2 
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+ 3- 
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(B.5) 



iV4 _ 2Ar3 - iV2 + 2Ar iV^ - GA^^ + llA^^ 



6A^ 
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Altogether, we obtain the following relations between the traces in the irreducible repre- 
sentations appearing in these decompositions and the traces in the fundamental representa- 
tion: 



HE 



I'^I'^IpI""! 
|33 



jj. 77] 



tr, 



\{{ii){y)^p){G) + {(/x)M(pa)} + m{ypG)}+ 

+ {(/XZ^)(p(7)} + {(/XZX/9(7)}) 

|((/.)H(p)(c7) + {(H(p)}(a) - {liG){y){p)- 
-{HG){up) - {fiG{u){p}) + {(/xi/p)}(cr) - {nG{up})) 

-(//)(i/Cr)(p) - {fMp){u){G) + {{fMu){pG)} - {npiy){G) - {fl){upG)- 
-{HGp){v) - {iXVg){p) - inGpu) - {pvpG) + {pGVp) + {ppVG)) 

\{{h){v){p){g) - {v)m{pG)} + {hu){p){g) - {pu){pG)- 

-{^iu{p){G}) + {p{pG}){l^) + {^iu{pG})) 

ir((M)M(p)M - m{uKpG)} + {ip){upG)}+ 

+ {{pu){pG)} - {{pupG)}) 



(B.6) 



Note that adding tr p for all the irreducible representations p with two and four boxes yields 



tr (r^)tr (r'^) 
tr (r'^) tr {r") tr (r^) tr {t") 



tr 
tr 
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(B.7) 



respectively, which confirms that our decomposition is correct. 

In the case in which all the factors of the Kronecker product are the same, the formulae 
given in ()B.6p reduce to character formulae, which are found e.g. in [31]. For convenience, we 
also list them in the following table: 



P 



m 

B 

rrm 



Xp(A) 



i((trA)2+ trA2) 
i((trA)2 - trA2) 



^ (( tr (A))"^ + 6( tr A)^ tr A^ + 8 tr A tr A^ + 3( tr A^)^ + 6 tr A"^) 
|(( tr (A))4 + 2( tr A)^ tr A^ - ( tr K^f - 2 tr A^) 
j^((tr (A))4 -4trAtrA3 + 3(trA2)2) 

|(( tr (A))4 - 2( tr A)^ tr A^ - ( tr K^f + 2 tr A^) 
^((tr (A))"^ - 6(trA)2trA2 + 8trAtrA3 + 3(trA2)2 - GtrA"^) 



C. Lie algebra conventions 

We use generators = {t^,t'^) of U(A^), which are orthonormal, i.e. 

tr (r^r'^) = 6^"" ; (C.9) 

= -^Iat, and are the Gell-Mann matrices generating SU(A^). In general, we use the 
Greek letters ^^v, p,a = 0, . . . , N"^ — 1 and Latin letters a, 6, c, d = 1, . . . , A^^ — 1 for labeling 
generators. The structure constants f°'^'^ and the symmetric tensor d"'^'^ are defined by the 
multiplication rule 

rV'' = ^,5"V° + ^iip^" + d^^'jr'^ . (C.IO) 
y N 2 



This implies that 



abc 



-itr ([r^rV^) 



jabc 



2 



N 



(C.ll) 



From the above definitions, it is obvious that f and d are totally antisymmetric and 
symmetric in their indices, respectively. Recall also the Fierz identity Tj^t^^ = Sudjk implying 
that r^r^ = iVlAr. 

The quadratic Casimir C2 with the standard eigenvalues l{l + 1) corresponds in our 
conventions to 

C2<I> = k[ra,[ra,n . (C.12) 

Starting from 



one easily proves the identity 



facefbde — J^i^a^cd — hc^ad) + {dbaedced — dbcedaed) 



(CIS) 
(C.14) 
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Furthermore, we have 

facedbde = ti {t^t^t't") - tr (rVVV'^) + tr (r W^r'') - tr (rVW) . (C.15) 
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